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It is currently a big challenge to accurately determine the symmetry energy Esymip) and the pure 
neutron matter equation of state -Epnm(p), even their values around saturation density po- We find 
that the electric dipole polarizability au in ^°®Pb can be determined uniquely by the magnitude of 
the Esymip) or almost equivalently the -Epnm(p) at subsaturation densities around po/3, shedding a 
light upon the genuine correlation between the q;d and the -Esym(p)- By analyzing the experimental 
data of the qd in ^°®Pb from RCNP using a number of non-relativistic and relativistic mean-held 
models, we obtain very stringent constraints on Esym{p) and -Bpnm(p) around po/3. The obtained 
constraints are found to be in good agreement with the results extracted in other analyses. In 
particular, our results provide for the hrst time the experimental constraints on i?pNM(p) around 
Po/3, which are in harmony with the recent determination of i?pNM(p) from microscopic theoretical 
studies and potentially useful in constraining the largely uncertain many-nucleon interactions in 
microscopic calculations of neutron matter. 
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1. Introduction .—The symmetry energy Esymip) 
well as the pure neutron matter equation of state (EOS) 
E'pnm(p) plays key roles in the investigation of physical 
objects from microscopic neutron-rich nuclei to macro¬ 
scopic neutron stars m and even in new physics beyond 
the standard model ■ Although significant progress has 
been made in recent years in understanding the Esymip) 
and Epnm(p) due to a lot of experimental, observa¬ 
tional and theoretical efforts, accurate determination of 
Esymip) and ifpNM(p), even their values around satura¬ 
tion density po ~ 0.16 fm“^, remains a big challenge (see, 
e.g.. Refs. [BHlI])- While heavy ion collisions and astro- 
physical observations provide two important approaches 
to constrain the symmetry energy from sub- to supra- 
saturation densities, nuclear structure probes usually can 
most effectively constrain the symmetry energy at sub¬ 
saturation densities. It has been established that nuclear 
mass can put stringent constraint on the magnitude of 
Esymip) around 2pol3 [T^l - fT7] and the neutron skin thick¬ 
ness Ar„p of heavy nuclei can fix the density slope L(p) 
of the symmetry energy around 2po/3 p^HTO] . At very 
low densities of O.OSpo < P < 0.2po and temperature in 
the range of 3 ^ 11 MeV where the clustering effects are 
essential, the symmetry energy has been obtained using 
data from heavy ion collisions pn| . 

In contrast to the Esymip), to our best knowledge, 
the only existing constraint on ifpNM(p) is EpnmIp = 
0.10 fm“^) = 11.4 ± 1.0 MeV obtained by analyzing the 
ground state properties of the doubly magic nuclei within 
the Skyrme-Hartree-Fock (SHF) approach [T^. Theo¬ 
retically, microscopic studies based on chiral effective 
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field theory (ChFFT) [21] [22] and Quantum Monte Carlo 
(QMC) calculations P5IB5] provide important informa¬ 
tion on ifpNM(p)- In these theoretical studies, the main 
uncertainty is due to the poorly known many-nucleon 
interactions. Therefore, experimental constraints on 
Fpnm(p) are extremely useful to understand the largely 
uncertain many-nucleon interactions in microscopic cal¬ 
culations of neutron matter. 

The nuclear electric dipole polarizability aD [26] has 
been proposed to be a good probe of the symmetry en¬ 
ergy m- However, their exact relationship has not yet 
been completely understood and even some controver¬ 
sial conclusions have been obtained in different analyses 
by examining the correlation between the ajo and the 
Esymip) around po (HHHj- In this work, we find that 
actually the au in ^°®Pb can be determined uniquely by 
the magnitude of the Esymip) or almost equivalently the 
Fpnm(p) at much lower densities around po/3, shedding 
a light upon the genuine correlation between the au and 
the Esymip)- This finding together with the q:d in ^°®Pb 
measured at the Research Center for Nuclear Physics 
(RCNP) [33] allows us to obtain quite precise constraints 
on Esymip) and i?pNM(p) around po/3. The present ex¬ 
perimental constraints on ifpNM(p) is potentially useful 
in constraining the poorly known many-nucleon inter¬ 
actions in the microscopic calculations of pure neutron 
matter. 

2. The symmetry energy and ajo ■—The EOS of asym¬ 
metric nuclear matter, defined by its nucleon specific 
binding energy, can be expanded as 

^ip, <J) = Eoip) -A Esymip)d^ + 0((5^), (1) 

where p = pn + Pp is nucleon density and <5 = (p„ — 
Pp)/iPp + Pn) is the isospin asymmetry with p„ (pp) de¬ 
noting the neutron (proton) density; Eq ip) represents the 
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EOS of symmetric nuclear matter; Egymip) is the sym¬ 
metry energy and it can be expressed as 


Esymip) 


1 d^E{p,6) 

2\ dP 


( 2 ) 


The Eq{p) is usually expanded around po as Eq[p) = 
Eo{po) + where the Kq is the 

so-called incompressibility coefficient. The E^ymip) can 
also be expanded around a reference density pr as 
E^y^[p) = E^yM + HPr){^) + 0((^)^) With 

L(pi.) = denoting the density slope of 

the symmetry energy at pr- Neglecting the higher- 
order terms in Eq. Q leads to the well-known empirical 
parabolic approximation ifpNM(p) ~ Eq{p) + Esy^ip)■ It 
should be emphasized that in the present work, all the 
results for Esymip) and ifpNM(p) are obtained exactly in 
mean-held models without parabolic approximation. 

The electric dipole polarizability od is proportional to 
the inverse energy-weighted sum of the electric dipole 
response [26] which is dominated by the isovector giant 
dipole resonance (IVGDR) — a nuclear collective oscil¬ 
lation of all the protons against all the neutrons with 
the symmetry energy E^ymip) acting as the restoring 
force [33]. Since the neutron and proton densities in the 
nuclear interior essentially do not change in the IVGDR, 
the od thus probes the symmetry energy not at po but 
rather at much lower densities around the nuclear surface 
where matter with extreme isospin or even pure neutron 
(proton) matter may form in the oscillation. 

A more quantitative preview about the an can be 
obtained from the macroscopic hydrodynamical model 
which predicts [35l |36] 


ao = 


24 



(3) 


where r represents the radial coordinate in the nuclei 
and Usym = E'syni(p)/p. Using an empirical 2-parameter 
Fermi distribution for the radial density distribution 
of nuclei and a simple parametrization of Esy^{p) = 
12.5(p/po)^/^-l-20(p/po)'’', one can find the in Eq. ^ 
is dominated by the symmetry energy values at low den¬ 
sities around po/3. Moreover, using a leptodermous ex¬ 
pansion in Eq. g, Lipparini and Stringari [35] derived a 
simple expression of od as 


ao 


12 




(4) 


where is the mean-square radius of the nucleus with 
mass number A, and by and bs are the volume and surface 
symmetry coefficients which are related to the symmetry 
energy coefficient asym(4l) of hnite nuclei with mass num¬ 
ber A as [33] 


^sym 


(^) 


1 by 

2 1 + (6«/6,)A-i/3- 
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Using Eq. Q and Eq. Roca-Maza et al. [331 recently 
obtained a relation between the od and the symmetry en¬ 
ergy parameters Esyya{po) and L{pq) by invoking the sim¬ 
ple relations by — ‘2Egyyy(^pQ^ and bsjby — ^E^yy^i^pQ'^jQ 
with Q being the surface stiffness coefficient [Tai33[3H]- 
On the other hand, substituting Eq. ([^ into Eq. Q, one 
can then obtain a new and interesting relation 


an (A) 


e" A(r") 

24 asym(]^A) 


( 6 ) 


which suggests that the od of a nucleus with mass num¬ 
ber A is inversely proportional to the symmetry energy 
coefficient of a nucleus with mass number (|)^A (e.g., 
for ^°®Pb, aD(A = 208) oc l/asym(A = 45) ). Consider¬ 
ing the strong correlation between the asym(A) and the 
EsymiPA) at a specific density pA [HI dll IMl SSjj one 
then expects that the an in ^°®Pb should be strongly 
correlated with Egymip) at p = pA =45 ~ Po/3 [40] . 

The above discussions indicate that a model- 
independent linear correlation may exist between I/od 
in ^*^®Pb and the magnitude of the symmetry energy 
around po/3. As will be shown in the following, this 
genuine correlation can be exactly confirmed by the mi¬ 
croscopic random-phase approximation (RPA) calcula¬ 
tions based on non-relativistic and relativistic mean-field 
models. Then given Esym(po/3) « Esym{po) - |T(po) or 
even better Asym(po/3) « Esyin(2po/3) - 4L(2po/3), one 
can easily understand the dependence of au in ^°®Pb on 
both Esymipo) and T(po) [IS] or on both ifsym(2po/3) 
and L(2po/3) [3^. This can also explain the relatively 
weak correlation of od in ^°®Pb with the Ar^p [33] since 
the latter is essentially determined by L(2po/3) [H] . 

3. The symmetry energy from od in ^^^Pb. —In or¬ 
der to study the correlation of the au in ^°®Pb with the 
symmetry energy value at different densities, we analyze 
the results from 62 representative non-relativistic and 
relativistic interactions which all give a good descrip¬ 
tion of the ground state properties of finite nuclei but 
predict very different density dependence of the symme¬ 
try energy, including 47 Skyrme interactions [131411143] 
(i.e., BSkl, BSk2, BSk6, Es, Gs, KDE, KDEOvl, MSk7, 
MSLO, MSLl, NRAPR, Rs, SAMi, SGI, SGII, SK255, 
SK272, SKa, SkI3, SkM, SkMP, SkM*, SkP, SkSl, SkS2, 
SkS3, SkS4, SkSGlS, SkT7, SkT8, SkT9, SKX, SKXce, 
SKXm, SkxslS, Skxs20, SLy4, SLy5, SLylO, v070, v075, 
v080, v090, vl05, vllO, Zs, Zs*) and 15 relativistic in¬ 
teractions involving FSU, NL3 and TF families [331 133] 
together with DD-ME2 [33] ■ For the Skyrme interactions, 
we evaluate the an using the Skyrme-RPA program by 
Colo et al. [33] while for the relativistic interactions, we 
directly invoke the results of RPA calculations reported 
in Refs. [SSISI]. 

The obtained data-to-data relations between IO^/od 
in ^°®Pb and the Esym{pr) at pr = 0.02, 0.05, 0.08, 0.11 
and 0.16 fm“^ are displayed in Fig. Also included in 
Fig.g are the linear fits together with the correspond¬ 
ing Pearson correlation coefficient r and the experimen- 
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FIG. 1: IO^/qd in vs. i?aym(pr) at pr = 0.02 fm“^(a), 0.05 fm“®(b), 0.08 fm“^(c), 0.11 fm“®(d) and 0.16 fni“®(e), 

predicted by a large number (62) of non-relativistic (SHF) and relativistic (RMF) models. The shaded regions correspond to 
linear fitting of IO^/qd = a + 6i5sym(Pr) with Icr and 2 (t confidence bands and r is the Pearson correlation coefficient. The red 
hatched band corresponds to the experimental result of q;d = 20.1 ± 0.6 fm® from RCNP |33 |. 


tal result of od = 20.1 ± 0.6 fm^ obtained from a high- 
resolution measurement at RCNP via polarized proton 
inelastic scattering at forward angles [33]. It is seen 
that the I/od exhibits a very strong linear correlation 
(r = 0.977) with the i?sym(p) around po/3 (i.e., pr = 0.05 
fm“^), confirming the predictions of both hydrodynam- 
ical and droplet models discussed earlier. The correla¬ 
tion remains strong at densities below po/S but it drops 
rapidly when the density is above about po/2. 

For Pr = 0.05 fm“^, the linear fit gives 

lOVttD = (18.94±0.85)-b(1.94±0.05)F;sym(pr), (7) 

with od in fm^ and £^sym(Pr) in MeV. Substituting the 
experimental value of od = 20.1 ± 0.6 fm^ into Eq. Q 





FIG. 2: Gonstraints on the symmetry energy iJaym(p) as a 
function of density p (see text for the details). The inset shows 
the density dependence of the Pearson correlation coefficient 
r between I/od in ^°®Pb and the Esymip)- 


then leads to 

EsymiPr) = 15.91 ± (0.77)exp ± (0.63)th, (8) 

where the uncertainties with “exp” and “th” are ob¬ 
tained from the propagation of the experimental un¬ 
certainty of ttD and the parameter errors in the lin¬ 
ear ht which reflect the systematic errors, respectively. 
Consequently, one can obtain a stringent constraint of 
E^ymiPr = 0.05 fm"^) = 15.91 ± 0.99 MeV. The simi¬ 
lar analysis allows us to extract constraints on the sym¬ 
metry energy value at other densities and the results 
from Pr = 0.02 fm“^ to 0.11 fm“^ are shown as red 
hatched band in Fig. [^ Shown in the inset in Fig. 
is the density dependence of the r value. It should be 
emphasized that in the present work, the symmetry en¬ 
ergy is obtained from the mean-field calculations without 
considering clustering effects. Indeed, at low densities 
(< Po/2) uniform nucleonic matter may become unsta¬ 
ble against cluster formation which may influence the 
symmetry energy. Theoretical studies [47j indicate that 
at zero temperature, only at very low densities (i.e., less 
than about 0.02 fm“^) where the fraction of light clus¬ 
ters becomes significant, clustering effects are essential 
and increase significantly the symmetry energy. There¬ 
fore, the constraints on the symmetry energy at densities 
below 0.02 fm“^ are not shown in Fig. [ 2 ] although the r 
value is still large. At higher densities (e.g., above 0.11 
fm“^), effective constraints cannot be obtained as the r 
value becomes much smaller. 

For comparison, we also shown in Fig.j^the constraints 
from transport model analyses of mid-peripheral heavy 
ion collisions of Sn isotopes (HIC) |48| and the SHF anal¬ 
yses of isobaric analogue states (IAS) as well as comb¬ 
ing additionally the neutron skin “data” (IAS-|-NSkin) 
in Ref. [15], and six constraints on the value of Esymip) 
around 2jipo from binding energy difference between 
heavy isotope pairs (Zhang) |16j . Fermi-energy difference 
in finite nuclei (Wang) [2], properties of doubly magic 
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nuclei (Brown) [T7], the giant dipole resonance in 
(Trippa) [49], the giant quadrupole resonance in 
(Roca-Maza) (SD) and the soft dipole excitation in ^^^Sn 
(Cao) |5T]. In addition, we also show the experimental 
results of the symmetry energies at densities below 0.2po 
and temperatures in the range 3 ~ 11 MeV from the anal¬ 
ysis of cluster formation in heavy ion collisions (Wada 
and Kowalski) |5D]. It is remarkable to see that a single 
data of od in ^°®Pb can give quite stringent constraints 
on Esym{p) around po/3 and the constraints are in very 
good agreement with other analyses. It is also interesting 
to see that the constrained Egymip) around po/7 is nicely 
consistent with the results extracted from heavy ion colli¬ 
sions (Wada) which consider the clustering effects [2D| . In 
addition, we note the i?sym(p) has been predicted in mi¬ 
croscopic calculations (see, e.g., Refs. [221[^IS3]) and the 
results are consistent with the experimental constraints 
shown in Fig. 

4- Neutron matter from od w —From the em¬ 

pirical parabolic approximation, one expects £'pnm(p) 
should play a similar role as £'sym(p) since i?o(p) is rel¬ 
atively well determined, especially around po/3. Using 
the similar analysis as in Fig. [^ we indeed find a strong 
correlation between l/ap in ^'^Pb and Upnm(p) around 
Po/3. Therefore, we can constrain £'pnm(p) around po/3 
in the similar way as constraining the i?sym(p) and the 
results are shown as red hatched band in Fig. [^ in the 
density interval 0.015 fm“^ < p < 0.11 fm“^ with the 
corresponding density dependence of r shown in the in¬ 
set. At lower densities (i.e., below 0.015 fm“^), although 
the r value remains close to unit and the clustering ef¬ 
fects are negligible in neutron matter, the constraints on 
A'pnm(p) are not shown in Fig.j^since the pairing effects, 
which are not considered in the mean-field calculations 
for Apnm(p), may become considerable [54]. We note 
the pairing effects on ifpNM(p) are negligibly small for 
p > 0.015 fm“^ [^. At higher densities (e.g., above 0.11 
fm“^), the r decreases rapidly and one cannot effectively 
constrain £'pnm(p)- 

Also shown in Fig. [^ are the predictions from ChEFT 
using N^LO potential in Ref. [3T] (ChEFT) and in 
Ref. [12] (Wellenhofer), Auxiliary-Field Diffusion Monte 
Carlo (AFDMC) calculations employing local N^LO 
ChEFT interaction with different cutoffs (QMC) [23], 
AFDMC calculations using adjusted nuclear force models 
by Gandolfi-Carlson-Reddy (GCR) [SS], Auxiliary-Field 
QMC calculations with the N^LO 2-body interactions 
plus the N^LO 3-body interactions (Wlazlowski) [24] . 
configuration interaction Monte Carlo calculations using 
nonlocal N^LO chiral interaction (Roggero) [IS], varia¬ 
tional calculations by Akmal-Pandharipande-Ravenhall 
(APR) [52], the Bethe-Bruckner-Goldstone calculations 
using many-body expansion up to the three hole-line 
level of approximation with quark model interactions 
with the auxiliary potential of gap choice (GC) or con¬ 
tinuous choice (BBG-QM 3h-gap and BBG-QM 3h- 
con) [SS], and the self-consistent Green’s function cal¬ 
culations that includes the effects of three-body forces 



p(fm-^) 


FIG. 3: Constraints and predictions for the pure neutron mat¬ 
ter EOS Epj^m{p) as a function of density p (see text for the 
details). The inset shows Pearson correlation coefficient r be¬ 
tween I/od in and i?pNM(p) as a function of p. 


(SCGF-N3LO-bN2LOdd) [^. At p = 0.1 fm-^, the con¬ 
straint ifpNM(p = 0.10 fm“^) = 11.4 ± 1.0 MeV is also 
shown (Brown) |17j . 

One can see from Fig. [^ that our present analyses of 
the data on au in ^°®Pb from RCNP give quite stringent 
constraints on £’pnm(p) around po/3. To our best knowl¬ 
edge, our present results provide for the first time the 
experimental constraints on E'pnm(p) around p « po/3. 
It is seen that our constraints are in excellent agree¬ 
ment with the predictions of APR, GCR and Wellenhofer 
as well as the constraint from Brown, and also consis¬ 
tent with other predictions. Interestingly, although our 
constraints are consistent with the predictions of both 
ChEFT and QMC within the uncertainty bands, there 
still exist some density regions where our constraints 
do not completely overlap with the uncertainty bands 
of ChEFT and QMC which are mainly due to the un¬ 
certainty of the many-body interactions. Therefore, our 
present experimental constraints on £'pnm(p) are poten¬ 
tially useful for constraining the many-body interactions 
in ChEFT and QMC calculations, which may be signifi¬ 
cant to improve the prediction at higher densities. 

5. Conclusion .—In summary, we have found that the 
electric dipole polarizability od in ^°®Pb can be deter¬ 
mined uniquely by the magnitude of the symmetry en¬ 
ergy Asyni(p) or almost equivalently the pure neutron 
matter EOS i?pNM(p) at subsaturaion densities around 
po/3, significantly deepening the understanding of the 
relation between the od and the Egymip)- This find¬ 
ing together with the an in ^°®Pb measured at RCNP 
has allowed us to obtain very stringent constraints on 
Asyin(p) and Apnm(p) around po/3. The present con¬ 
straints should be less model dependent since they are 
based on a large set of both non-relativistic and rela¬ 
tivistic models. Our results provide for the first time the 
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experimental constraints on -E'pnm(p) around po/3 which 
are potentially useful in constraining the many-nucleon 
interactions in microscopic calculations of neutron mat¬ 
ter. 
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